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A  NONCONFORMING  APPROXIMATE  SOLUTION  TO  A 
SPECIALLY  ORTHOTROPIC  AXISYMMETRIC 
THIN  SHELL  SUBJECTED  TO  A  HARMONIC 
DISPLACEMENT  BOUNDARY  CONDITION 

1.  INTRODUCTION 

The  design  of  lightweight,  inexpensive,  and  strong  shells  is  critical  to  many 
mechanical  systems.  The  analysis  of  various  types  of  shells  under  different  loading 
configurations  has  been  reported  by  numerous  researchers  (Naghdi  and  Berry,  1954; 
Flugge,  1960;  Junger  and  Feit,  1986;  Vinson  and  Sierakowski,  1987;  Qaisi,  1989).  In 
these  papers,  the  governing  partial  differential  equations  to  evaluate  a  particular  shell  are 
usually  developed  and  then  a  response  is  formulated.  The  response  is  an  equation  of 
motion,  which  is  required  to  compare  test  data  to  analytical  models  or  to  solve  an  inverse 
problem  (i.e.,  the  test  is  run  and  the  system  parameters  are  identified  by  the  test  data.) 

In  this  report,  the  partial  differential  equations  of  an  axisymmetric  composite  thin 
shell  subjected  to  a  harmonic  boundary  condition  are  developed.  First,  the  model  of  a  shell 
subjected  to  axial  tension  with  a  mechanical  shaker  at  the  forward  end  and  terminated  by  a 
rope  to  ground  at  the  after  end  is  formulated.  The  partial  differential  equations  are 
transformed  into  ordinary  differential  equations.  An  approximate  solution  to  the  ordinary 
differential  equations  of  motion  is  then  derived.  The  model  is  verified  with  comparisons  to 
results  from  finite  element  analysis  of  the  same  thin  shell  problem.  The  model  is  also 
compared  to  a  previous  model  of  the  system. 

The  new  model  corresponds  to  the  testing  configuration  in  the  Axial  Vibration  Test 
Facility  (AVTF)  at  the  Naval  Undersea  Warfare  Center  (NUWC),  as  shown  in  figure  1. 
This  facility  has  been  designed  to  provide  a  simple  procedure  for  testing  shells  under 
varying  tensions  and  temperatures. 
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Figure  1.  Laboratory  Configuration 


2.  SYSTEM  MODEL 

The  system  model  is  an  axisymmetric  thin  shell  with  a  displacement-driven 
boundary  condition  at  one  end  and  a  mechanically  grounded  spring  at  the  other  end.  The 
governing  partial  differential  equations  are  modified  axisymmetric  Love-Timoshenko 
formulas.  The  longitudinal  (axial)  partial  differential  equation  that  relates  the  shells 
longitudinal  acceleration  to  its  axial  stress  is 

_d^u{x,t)  dar{x,t) 

where  u{x,t)  is  the  axial  displacement  (m).  x  is  the  axial  spatial  location  (m),  t  is  time  (s),  p 
is  the  density  of  the  shell  (kg/m^),  and  Oxix.t)  is  the  axial  stress  in  the  shell  (N/m^).  The 
radial  (circumferential)  partial  differential  equation  that  relates  the  radial  acceleration  of  the 
shell  to  its  tension,  bending  stiffness,  and  circumferential  stress  is 


d^w{x,t)  T  d^w{x,t)  „  d'^wjxj)  h 


2na  dx^ 


—  aAxj)  , 


where  w{x,t)  is  the  radial  displacement  (m),  h  is  the  shell  thickness  (m),  T  is  the  tension  in 


the  shell  (N),  a  is  the  average  radius  of  the  shell  (m),  B  is  the  bending  stiffness  of  the  shell 
(Nm),  and  cr^x.r)  is  the  circumferential  (hoop)  stress  (N/m^).  The  bending  stiffness  of 


the  shell  is  given  by  the  expression 


B  = 


(3) 
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where  Ex  is  the  axial  elastic  modulus  (N/m^).  The  longitudinal  elastic  modulus  of  a  thin 
shell  is  related  to  Young's  modulus  using  the  equation 

— - - T.  (4) 

where  E  is  Young's  modulus  (N/m^),  Ox0is  Poisson's  ratio  relating  an  axial  force  to  a 
circumferential  displacement,  and  is  Poisson's  ratio  relating  a  circumferential  force  to 
an  axial  displacement.  The  relationship  between  these  two  Poissons'  ratios  will  be 
discussed  later  in  the  report. 

The  dynamic  stresses,  <Tx{x,t)  and  <T0(x./),  are  determined  by  combining  two- 
dimensional  stress-strain  relations  and  the  strain-displacement  relations  in  cylindrical 
coordinates.  The  equation  relating  the  longitudinal  stress  to  the  axial  and  circumferential 
strain  for  an  axisym metrical  solid  is  (Jones,  1975) 


CTjt(x,r)  =  Ej^£x{x,t)+  Vx(^E^e^{x,t)  ,  (5) 

where  e;cU,0  is  the  longitudinal  strain  (dimensionless),  E(p  is  the  circumferential  modulus 
(N/m2),  and  f^x.r)  is  the  circumferential  strain  (dimensionless).  The  longitudinal  strain  is 
related  to  the  longitudinal  displacement  by 
du{x,t) 


£xix,t)  = 


dx 


(6) 


and  the  axisymmetric  circumferential  strain  is  related  to  the  circumferential  displacement  in 
a  thin  shell  by 

w(x,f) 


£<l>{x,t)  = 


(7) 


Equations  (6)  and  (7)  can  be  inserted  into  equation  (5),  which  produces  an  axial  stress  of 


r.  Mx,t)  .  Vx<l>E0W{x,t) 

iKX,t)-tx — - - + -  . 

ax  a 


(8) 


Tlie  equation  relating  the  hoop  stress  to  the  circumferential  and  longitudinal  strain  for  an 
axisymmetrical  solid  is  (Jones,  1975) 
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a^ix,t)  =  E0e^ix,t)+  .  (9) 

Equations  (6)  and  (7)  can  be  inserted  into  equation  (9)  to  yield  a  circumferential  stress  of 
E^w(x,t)  duixj) 

(T^ix,t)  =  -^ - +  .  (10) 

Equation  (8)  is  now  differentiated  with  respect  to  x  and  inserted  into  equation  (1), 


which  gives 


d^u{x,t)  d^uixj)  ^x<l)^0  dw{x,t) 


Equation  (10)  is  next  inserted  direcUy  into  equation  (2)  yielding 

d^w{x,t)  ^d'^wixj)  ,  T  d^w{xj)  hE^  du(,x,t) 

ph — —  =  -fi — 71—  +  :; - 71 - f^w(jr,r) - - — .  (12) 

dr  dx^  2na  dx  o  a  dx 

Equations  (11)  and  (12)  are  the  two-dimensional  axisymmetric  equations  of  motion  for  the 

specially  orthotropic  shell. 

The  harmonic  boundary  condition  in  the  axial  direction  at  jr  =  0  is  modeled  as 

«(0,0  =  (/o^'‘"  ,  (13) 

where  Vq  is  the  amplitude  of  the  displacement  at  the  boundary  (m),  i  is  the  square  root  of 
-1,  and  (o  is  frequency  (rad/s).  The  boundary  condition  in  the  axial  direction  at  x  =  L  is 
modeled  as  a  spring  connected  at  one  end  to  the  shell  and  to  mechanical  ground  (zero 
displacement)  at  the  other  end.  This  equation  is  derived  as  follows  by  equating  the  force  at 
the  end  of  the  shell  to  the  force  in  the  spring: 

=  (14) 

dx 

where  A  is  the  cross-sectional  area  of  the  shell  (m^)  and  ^  is  the  spring  constant  (N/m). 
The  boundary  conditions  in  the  radial  direction  are  zero  displacement  and  zero  shear  force 
(lunger  and  Feit,  1986).  These  conditions  are  written  as 

MO, 0  =  0,  (15) 

w(L,0  =  0,  (16) 
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and 


dx 

dw{L,t) 

dx 


=  0  . 


=  0  . 


(17) 


(18) 


Although  the  boundary  conditions  given  by  equations  (15)  -  (18)  are  not  used  in  the  final 
analysis,  they  are  included  here  for  a  complete  treatment  of  the  theory.  The  effect  of  not 
using  these  boundary  conditions  is  discussed  later  in  the  report. 

It  should  be  noted  that  the  previous  equations  used  to  model  this  system  are  (Chase, 

1975) 


d^uixj)  _  d^u{x,t)  vaEr  dw{x,t) 


dt^ 


=  E 


a 


dx 


(19) 


and 


d^w{x,t)  ^d^wix,t)  T  d^wix,t)  ha^Ex  ,  ,  hvaEx  duix,t) 

p*- Ti- = - :5-  - 


dt^ 


dx 


a 


dx 


(20) 


with 


2 

a  =  — 


(21) 


where  a  and  are  both  dimensionless  quantities.  Equations  (11)  and  (12),  derived 
earlier  in  this  section,  will  be  compared  to  equations  (19)  and  (20)  as  well  as  to  results  from 
finite  element  analysis  of  several  thin  shell  problems. 


3.  A  NONCONFORMING  APPROXIMATE  SOLUTION 
An  approximate  closed  form  solution  can  be  developed  for  the  partial  differential 
equations  and  the  boundary  conditions  in  the  longitudinal  direction.  This  calculation 
assumes  that  the  time  and  space  modes  of  the  system  will  decouple.  A  steady  state  solution 
is  also  of  interest,  because  it  allows  the  displacements  to  be  written  as 

uix,t)  =  U{x)e'^  (22) 
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and 

w{x,t)=W{x)e‘‘^  .  (23) 

A  second  assumption,  which  presumes  that  the  dynamic  coupling  from  the  radial  equation 
into  the  longitudinal  equation  is  negligible,  results  in  an  approximate  solution  to  the 
problem.  It  will  be  shown  that,  based  on  the  selected  loading  condition  on  the  structure, 
this  premise  is  reasonable.  If  the  radial  displacement  derivative  term  is  dropped  and 
equation  (22)  is  used,  equation  (11)  can  be  rewritten  as 


with 


d^U{x)  ,2,,/  s  rx 
dx 


(24) 


03 


(25) 


where  k  is  the  extensional  wavenumber  (rad/m).  The  general  solution  to  equation  (24)  is 
(/(x)  =  Geos itr  +  Z/ sin  fee  .  (26) 

The  boundary  conditions  from  equations  (13)  and  (14)  are  now  applied  to  equation  (26)  to 
solve  for  the  constants  G  and  H  resulting  in 


Uix) 

Uo 


=  cos  foe + 


AE^ksinkL  -  ^  coskL] 
AEjrk  cos  kL  +  ^  sin  kL  J 


sin  foe  , 


(27) 


which  is  the  closed  form  approximate  solution  to  the  axial  equation  of  motion. 

The  radial  equation  of  motion  (12)  is  now  rewritten  using  equations  (22),  (23),  and 
the  first  derivative  of  equation  (27)  as 


d'^Wix) 

T 

d'^Wix)  ( hE^ 

dx* 

llta 

H'(r)  = 


- - - (sin fox- //cos fox)  , 


(28) 


where 
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, ,  AE^k sin  kL  -  ^  cos  kL 

H  =  -  . 

AEjrk  coskL  +  q  sin  kL 


(29) 


Equation  (28)  is  now  a  nonhomogeneous  ordinary  differential  equation  whose  forcing 
function  was  derived  from  the  axial  equation  of  motion.  The  particular  solution  to  equation 
(28)  is  of  the  form 

W(x)  =  Mcoskx  + Nsinkx  .  (30) 

Equation  (30)  can  be  differentiated  two  and  four  times  with  respect  to  x  and  then  inserted 
into  equation  (28),  yielding  the  radial  displacement  of 


W(x)  (Inahv^j^Ej^k)  (sinitx  -  Hcoskx) 

Uq  ~  {Ina^Bk'^  +  Tak^  -^InhE^-lna^co^ph) 


(31) 


The  boundary  conditions  on  the  radial  equations  are  not  used  to  define  the  equations  of 
motion  in  the  radial  direction,  therefore  the  circumferential  solution  does  not  conform  to  the 
boundary  conditions  at  j:  =  0  and  x  =  L.  Others  (Vinson  and  Sierakowski,  198?)  note  that 
the  homogeneous  solution  to  shell  equations  of  motion  goes  to  zero  away  from  either  end 
of  the  shell.  Equations  (27)  and  (31)  now  define  the  displacements  of  the  shell  in  the  axial 
and  radial  directions. 


4.  FINITE  ELEMENT  ANALYSIS 

The  closed  form  solutions  (equations  (27)  and  (31))  are  compared  to  results 
obtained  using  finite  element  analysis.  Finite  element  analysis  is  a  discretized  modeling 
technique  that  breaks  down  the  structure  into  a  number  of  subdomains  (elements)  (Cook, 
1974;  Zienkiewicz,  1983).  Constitutive  equations  are  applied  to  each  of  these  elements  to 
develop  element  equations  of  motion.  These  individual  equations  of  motion  are  then 
assembled  to  form  a  global  equation  of  motion.  The  global  equation  of  motion  is  of  the 
form 


M 


dt^ 


dt 


+  Kg,(r)  =  f(x,r)  , 


(32) 
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where  M  is  the  system  mass  matrix,  C  is  the  system  damping  matrix,  K  is  the  system 
stiffness  matrix,  g,Y t)  is  the  generalized  displacement  vector  in  the  /th  direction,  and  f  is  the 
external  forcing  function  vector.  The  finite  element  results  displayed  here  were  obtained 
using  ANSYS,  which  is  a  general  purpose  finite  element  program.  Finite  element  analysis 
has  the  advantage  of  handling  material  and  property  discontinuities  that  are  not  easily 
resolved  with  other  numerical  techniques.  It  is  also  used  for  comparisons  with  other 
experiments  and  theories  to  verify  results.  ANSYS  has  shell  models  that  are  internal; 
therefore,  the  equations  of  motion  derived  here  (equations  (11)  and  (12))  were  not  input 
into  the  finite  element  program  as  a  basis  for  analysis  (Kohnke,  1987).  Although  neglected 
in  the  development  of  equations  (27)  and  (31),  the  effect  of  dynamic  coupling  from  the 
radial  equation  into  the  longitudinal  equation  is  included  in  the  finite  element  analysis.  It  is 
important  to  note  that  the  finite  element  analysis  is  an  approximate  numerical  technique;  it  is 
extremely  different  from  the  analytical  solution  presented  here.  A  major  disadvantage  of 
finite  element  analysis  is  its  inability  to  yield  clo.sed  form  solutions;  thus,  any  time  a 
parameter  is  changed,  the  entire  problem  must  be  resubmitted  for  analysis.  It  is  also 
numerically  intensive  and  cannot  be  implemented  easily  on  small  computers. 

5,  AN  ISOTROPIC  EXAMPLE 

The  closed  form  solutions  derived  earlier  are  now  analyzed  using  an  isotropic 
example.  Isotropic  shells  have  the  elastic  modulus  and  Pois.son’s  ratio  in  the  axial  direction 
equal  to  the  circumferential  direction.  The  shell  properties  are  listed  in  table  1.  These 
closed  form  solutions  are  compared  to  results  derived  using  finite  element  analysis. 
Equations  (19)  and  (20)  are  identical  to  equations  (11)  and  (12)  for  the  isotropic  example 
and  thus  are  not  considered  separately  for  this  case.  Figure  2  shows  the  model  results 
compared  to  the  finite  element  analysis  for  the  longitudinal  ((/)  displacement,  and  figure  3 
depicts  the  model  results  compared  to  finite  element  analysis  for  the  radial  (VP) 
displacement.  The  closed  form  solutions  for  both  the  radial  and  longitudinal  displacements 
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are  in  very  good  agreement  with  the  finite  element  analysis.  Also,  as  predicted  by 
equations  (27)  and  (31),  the  radial  resonances  are  driven  by  the  longitudinal  resonances. 


Table  1.  Shell  Properties  for  IsoU'opic  Example 


Geometric  Properties 


J 


Length,  L 

25  m  (82  ft) 

Respon.se  location,  x 

7  m  (23  ft) 

Shell  radius,  a 

0.0254  m  ( 1  in) 

Shell  thickness,  h 

0.(K)254  m  (0.1  in) 

Shell  cross-sectional  area,  A 

0.0(K)4054  m2  (0.6284  in2) 

Material 

Properties 

Longitudinal  modulus.  Ex 

1  X  109  N/m2 

Circumferential  modulus,  £0 

1  X  109  N/m2 

Poisson's  ratio,  Vx0 

0.49  (dimensionless) 

Poisson's  ratio,  V0x 

0.49  (dimensionless) 

Density,  p 

1000  kg/m-^ 

Testing  Properties 

Spring  constant  ^  at  x  =  L 

1  X  105  N/m 

Shell  tension,  T 

4450  N (1000  lb) 

9 


Frequency  (Hz) 


Figure  3.  Radial  Response  of  Isotropic  Example  at  x  =  7  Meiers 
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FINITE  LENGTH  SHELL 


TR  10,165 


The  closed  form  solutions  are  now  compared  to  a  specially  orthotropic  example. 
Specially  orthotropic  describes  the  property  where  the  elastic  modulus  and  Poisson's  ratio 
are  different  in  the  radial  direction  than  in  the  axial  direction.  This  property  also  requires 
that  the  major  axis  of  the  material  properties  he  aligned  with  the  circumferential  and 
longitudinal  directions  of  the  shell.  Specially  orthotropic  shell  properties  are  listed  in  table 
2.  The  material  properties  of  a  specially  orthotropic  composite  satisfy  the  equation  (Jones, 
1975) 


V 


_  ^0x 
’■“x 


(33) 


From  equation  (33),  it  is  evident  that  the  material  possesses  only  three  independent 
parameters.  For  comparison  to  equations  (27)  and  (31),  the  approximate  solution 
technique  developed  in  section  3  can  be  applied  to  equations  (19)  and  (20).  The  solution  in 
the  axial  direction,  which  is  the  same  as  that  developed  earlier,  is  given  by  equation  (27). 
The  equation  of  motion  in  the  radial  direction  is 

lY(x)  _  ilnahvaEjk)  {sinkx-  Hcoskx) 


Uo 


(Irca^Bk^  +  Tak^  +  InhE^  -  2na^O)"ph) 


(34) 


where  equation  (34)  is  the  solution  to  equation  (20).  Equation  (34),  although  similar  to 
equation  (31),  results  in  a  much  different  radial  response  of  the  structure  for  a  specially 
orthotropic  example.  Figure  4  shows  the  model  results  (equation  (27))  compared  to  finite 
element  analysis  for  the  longitudinal  ((/)  displacement,  and  figure  5  illustrates  the  model 
results  (equation  (31))  compared  to  finite  element  analysis  and  the  model  given  in  equation 
(34)  for  the  circumferential  (VV)  displacement  with  v  =  used  in  equation  (34).  Figure  6 

presents  the  model  results  (equation  (31))  compared  to  finite  element  analysis  and  the 
model  given  in  equation  (34)  for  the  circumferential  (vt')  displacement  with  i>  =  V0x  in 

equation  (34).  The  comparison  for  equation  (34)  uses  both  values  of  Poisson's  ratio 
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because  the  original  reference  (Chase,  1975)  is  unclear.  It  is  shown  that  equation  (34) 
either  overestimates  or  underestimates  the  structural  response  depending  on  the  value  of 
Poisson's  ratio  used.  The  phase  angle  of  equation  (34)  is  not  shown  in  figures  5  and  6;  it 
is  the  same  as  the  phase  angle  of  equation  (31),  therefore  it  is  omitted. 


Table  2.  Shell  Properties  for  Specially  Orthotropic  Example 


Geometric  Properties 

Length,  L 

25  m  (82  ft) 

Response  location,  x 

16  m  (52.5  ft) 

Shell  radius,  a 

0.0254  m(l  in) 

Shell  thickness,  h 

0.(K)254  m  (0.1  in) 

Shell  cross-sectional  area,  A 

0.0004054  m2  (0.6284  in2) 

Material 

Properties 

Longitudinal  modulus.  Ex 

1  X  109  N/m2 

Circumferential  modulus, 

1  X  1()2  N/m2 

Poisson’s  ratio,  Vx0 

0.4900  (dimensionless) 

Poisson's  ratio,  V0x 

0.(X)49  (dimensionless) 

Density,  p 

1000  kg/m-' 

Testing 

Properties 

Spring  constant  ^  at  x  =  L 

1  X  105  N/m 

Shell  tension,  T 

2225  N (500  lb) 
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Frequency  (Hz) 


Figure  4.  Longitudinal  Response  of  Specially  Orthotropic  Example  atar  =  16  Meters 
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Phase  of 


X 

^  Finite 
Element 
Solution 


Closed  Form 
Solution 
Equation  (31) 


75  1C 

Frequency  (Hz) 


. I  ......  I  I  .....  .  .  , . , 

I  25  50  75  100 

Frequency  (Hz) 

Figure  6.  Radial  Response  of  Specially  Orthotropic  Example  at 
X  =  16  Meters  Using  v  =  0.4900  for  Equation  (34) 
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7.  A  SPECIALLY  ORTHOTROPIC  EXAMPLE  - 
SEMI-INFINITE  LENGTH  SHELL 

A  semi-infinite  specially  orihotropic  shell  is  now  analyzed.  This  type  of  analysis  is 
important  when  the  shell  is  extremely  long  or  when  all  the  energy  moving  down  the  shell  is 
absorbed  either  internally  or  at  the  boundary.  A  convenient  way  to  alter  equations  (1) 
through  (18)  into  a  semi-infinite  system  is  to  change  the  boundary  condition  dl  x  -  L 


(equation  (14))  to 


3u(L,t)  du(L,t) 


where  c  now  corresponds  to  a  viscous  damper  (Ns/m)  at  the  end  of  the  shell  as  depicted  in 
figure  7.  It  can  be  shown  that  a  value  of  (Hull.  1993) 

c  =  A^[^  (36) 

produces  a  longitudinal  response  that  contains  only  propagating  wave  energy.  Although 
the  shell  has  finite  length,  all  the  energy  moving  down  the  shell  will  be  absorbed  by  the 
viscous  damper,  and  the  shell  response  appears  as  if  it  is  infinitely  long.  Inserting  equation 
(36)  into  equation  (35)  and  solving  the  equations  of  motion  ((11)  and  (12))  with  the 
boundary  condition  of  equation  (13)  yields  the  longitudinal  response  of 

(37) 

Vo 

and  the  radial  response  of 

Wjx)  ^ _ {27tahv^^E^ik)e-'^ 

Uq  ~  {Ina^Bk'^  +  Tak^  2nhE^  -  Ina^co^ph)  ' 

Equations  (19)  and  (20)  (Chase,  1975)  can  also  be  solved  using  the  same  viscous 
boundary  at  jc  =  L.  The  axial  equation  of  motion  is  identical  to  equation  (37).  The  radial 
equation  of  motion  is 

Wix)  {2nahvaE^ik)e-'^ 

Uq  {2na^Bk^  Tak^ +2nhE^  -  2na^o}^ph) 
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Impedance 

i  vi’(x.r) 

Impedance 

Head 

I 

Head 

1 

Shaker  ■ 

a 

— y. 

- ►x.  u(x.r) 

Shell 

Damper 

c 

Figure  7.  Configuration  of  Semi-Infinite  System 


The  closed  form  solutions  are  now  compared  for  the  semi-infinite  specially 
orthotropic  shell.  The  shell  properties  are  listed  in  table  7>.  Figure  8  shows  the  model 
results  (equation  (37))  compared  to  finite  element  analysis  results  for  the  longitudinal  (U) 
displacement.  Although  a  finite  length  is  used  in  the  finite  element  analysis,  the  value  is 
also  arbitrary.  Figure  9  presents  the  model  re.sults  (equation  (38))  compared  to  finite 
element  analysis  results  and  the  model  given  in  equation  (39)  for  the  circumferential  (VV) 
displacement  with  v  =  Vx0.  Figure  10  are  the  model  results  (equation  (38))  compared  to 

finite  element  analysis  results  and  the  model  given  in  equation  (39)  for  the  circumferential 
(W)  displacement  with  v  =  v^x-  The  closed  form  solutions  given  by  equations  (37)  and 

(38)  agree  very  well  with  results  from  finite  element  analysis.  As  in  the  previous  example, 
equation  (39)  either  overestimates  the  response  or  underestimates  the  response,  depending 
on  the  value  of  Poisson’s  ratio  used. 
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Table  3.  Shell  Properties  for  Specially  Onhotropic  Semi-Infinite  Length  Example 


— 

n  Geometric  Properties 

Length,  L 

Semi-infinite 

Response  location,  x 

12  m 

Shell  radius,  a 

0.0254  m  ( 1  in) 

Shell  thickness,  h 

0.00254  m  (0.1  in) 

Shell  cross-sectional  area,  A 

0.0004054  m2  (0.6284  in2) 

Material 

Properties 

Longitudinal  modulus.  Ex 

1  X  11)9  N/m2 

Circumferential  modulus,  E(p 

1  X  107  N/m2 

Poisson’s  ratio, 

0.4900  (dimensionless) 

Poisson's  ratio, 

0.(K)49  (dimensionless) 

Density,  p 

1000  kg/m^ 

Testing  Properties 

Damper  constant  c  at  x  =  L 

405.4  Ns/m 

Shell  tension,  T 

2225  N  (500  lb) 
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Example  at  jc  =  12  Meters 
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Figure  9,  Radial  Response  of  Semi-Infinite  Specially  Orthotropic  Example  at 
x=  \2  Meters  Using  v  =  0.0049  for  Equation  (39) 
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Figure  10.  Radial  Response  of  Semi*Infinite  Specially  Orthotropic  Example  at 


x=  12  Meters  Using  v  =  0.4900  for  Equation  (39) 
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8.  A  SPECIALLY  ORTHOTROPIC  EXAMPLE  - 
INCLUSION  OF  STRUCTURAL  DAMPING 

The  model  can  include  structural  damping  by  changing  the  elastic  modulus  from  a 
real-valued  term  to  a  complex-valued  term.  The  damping  is  modeled  as  an  (imaginary)  loss 
modulus  with  the  (real)  storage  modulus.  The  shell  properties  for  this  damped  example  are 
listed  in  table  4.  Except  for  the  addition  of  damping,  these  properties  are  the  same  as  those 
used  for  the  example  in  section  6.  Figure  1 1  shows  the  damped  model  results  compared  to 
the  undamped  model  results  for  the  longitudinal  (U)  displacement,  and  figure  12  illustrates 
the  damped  model  results  compared  to  the  undamped  model  results  for  the  circumferential 
(WO  displacement.  ANSYS  does  not  have  the  capability  to  model  specially  orthotropic 
materials  with  varying  damping  coefficients  in  different  directions;  therefore,  no 
compari.son  to  finite  element  analysis  is  included  here.  If  equation  (20)  were  compared  to 
this  example,  the  difference  would  be  similar  to  that  shown  in  figures  5  and  6.  The  results 
from  figures  1 1  and  12  show  that  the  addition  of  damping  reduces  the  levels  of  shell 
resonance  and  increases  the  levels  of  shell  antiresonance.  The  transfer  function  is 
otherwise  very  similar  to  the  undamped  example.  The  effect  of  damping  would  be  greater 
if  the  shell  were  longer,  because  the  point  harmonic  forcing  function  would  be  absorbed  by 
a  larger  medium. 
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Table  4.  Shell  Properties  for  Specially  Onhotropic  Damped  Example 


Geometric  Properties 

Length,  L 

25  m 

Response  location,  x 

16  m 

Shell  radius,  a 

0.0254  m  (1  in) 

Shell  thickness,  h 

0.00254  m  (0.1  in) 

Shell  cross-sectional  area,  A 

0.0004054  m2  (0.6284  in2) 

Material  Properties 

Longitudinal  modulus.  Ex 

1  X  109(1-h().05/)  N/m2 

Circumferential  modulus, 

1  X  107(l-i4).20/)N/m2 

Poisson's  ratio,  Vx((t 

().49(K)  (dimensionless) 

Poisson's  ratio,  v^x 

0.0049  (dimensionless) 

Density,  p 

KXK)  kg/m3 

Testing  Properties 

Spring  constant  ^  at  x  =  L 

1  X  105  N/m 

Shell  tension,  T 

2225  N (500  lb) 
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Frequency  (Hz) 

Figure  1 1.  Longitudinal  Response  of  a  Damped  and  Undamped  Specially 
Onhotropic  Example  at  jc  =  16  Meters 
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9.  DISCUSSION 


All  the  shells  modeled  in  this  study  compare  favorably  to  finite  element  theory.  The 
approximate  closed  form  solution  is  very  accurate  for  a  harmonic  displacement  boundary 
condition.  It  is  possible  that  this  technique  will  not  be  as  accurate  for  other  types  of  forcing 
or  other  geometric  sizes.  The  equations  previously  used  to  model  this  system  (equations 
(19)  and  (20),  Chase,  1975)  do  not  produce  valid  results  for  specially  orthotropic  analysis. 
They  can,  however,  be  easily  reformulated  to  conform  to  the  analysis  given  here.  The 
radial  coupling  term  in  equation  (19)  needs  to  be  altered  to 


2 

VCcEx  dw(x,t)  ^ 
a  dx  a  dx 

and  the  longitudinal  coupling  term  in  equation  (20)  to 


(39) 


hvaEjr  du{x,t)  du{x,t)  . 

a  ax  a  ax 

Such  a  modification  will  produce  results  that  are  in  agreement  with  the  finite  element 
analysis.  Note  that  the  larger  the  difference  between  the  axial  and  circumferential  modulus, 
the  greater  is  the  discrepancy  between  equations  (12)  and  (20).  When  the  axial  modulus 
equals  the  circumferential  modulus,  equation  (12)  is  equal  to  equation  (20). 

The  nonconforming  solution  was  compared  to  the  finite  element  solution  near  the 
ends  of  the  shell  at  x  =  0.5  m  and  x  =  24.5  m  for  the  examples  listed  in  sections  5,  6,  and 
7.  These  are  the  two  locations  in  the  finite  element  model  that  are  closest  to  the  end  of  the 
shell  where  a  solution  can  be  obtained.  The  agreement  between  the  finite  element  model 
and  the  closed  form  approximate  solution  was  similar  to  that  shown  earlier.  This 
comparison  shows  that  the  radial  boundary  conditions  do  not  influence  the  displacements 
of  the  shell  even  0.5  meters  near  its  end.  If  the  nonconforming  solutions  were  compared  to 
the  finite  element  results  at  x  =  0  and  x  =  L,  the  solutions  would  be  similar  in  the 
longitudinal  direction  and  dissimilar  in  the  radial  directions.  This  dissimilarity  would  only 
be  a  concern  if  the  solution  very  close  to  the  end  of  the  shell  was  of  interest 
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10.  CONCLUSIONS 

A  nonconforming  approximate  solution  to  a  thin  shell  subjected  to  a  harmonic 
displacement  boundary  condition  was  developed  in  this  report.  This  closed  form  solution 
agrees  well  with  finite  element  solutions  for  the  several  models  evaluated.  It  is  shown  that 
previous  governing  equations  for  the  case  of  specially  orthotropic  shell  models  are 
erroneous.  These  equations  can  be  reformulated  to  match  the  analysis  here  by  an  alteration 
in  the  coupling  terms  between  the  radial  and  longitudinal  equations  of  motion.  Future  work 
in  this  area  should  include  laboratory  testing  to  obtain  experimental  results. 
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